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Wc study B^ J/tpTT^ and B^ J/ipK^ decays in a sample of about 89 million BB pairs 
collected with the BABAR detector at the PEP-II asymmetric B-factory at SLAC. We observe a signal 
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of 244 ± 20 ^ J/i)-K^ events and determine the ratio B{B^ J/ipTv^)/B{B^ J/ipK^) to be 
[5.37±0.45(stat.)±0.11(syst.)l%. The charge asymmetries for the ~* J/%1)-k^ and B^ ~* J/ip 
decays are determined to be An = 0.123±0.085(stat.)±0.004(syst.) and Ak = 0.030±0.015(stat.) ± 
0.006(syst.), respectively. 

PACS numbers: 13.25.Hw, 12.38.Qk 



We present an analysis of J/ipir^ and 

J/ipK^ decays that measures the ratio of branching frac- 
tions and searches for direct CP violation. The Cabibbo- 
suppressed decay B^ J/ipn^ proceeds via a, b ^ qd 
transition. It is expected to have a rate about 5% of 
that of the Cabibbo-allowed mode B^ J/%l;K^. The 
Standard Model predicts that for b ^qs decays the tree 
and penguin contributions have the same weak phase and 
thus no direct CP violation is expected in B^ J/^ 
decays. However, for b — *■ gd, the tree and penguin 
contributions have different phases and charge asym- 
metries as large as a few percent may occur 0. In 
the absence of isospin violation, the CP asymmetry in 
B± ^ J/tJjK^ provides a measurement of the ratio 
where A (A) is the decay amplitude for the neutral 
mode B"{B°) J/^K^. 

Previous studies of the B^ J/tjj tt^ mode have been 
performed by the CLEO 0, CDF 0, BABArB and 
Belle la] collaborations. The PDG 2002 average 01 of the 
ratio of branching fractions is (4.2±0.7)%. A recent Belle 
result gives ^ J/V'7r±) = (3.8 ± 0.6 ± 0.3) x 10"^. 

The PDG 2002 averages of the charge asymmetries are 
= -0.01 ±0.13 and Ak = -0.007 ±0.019 (see Eq.H 
for the definition of the sign of the asymmetry) . 

The analysis reported in this paper is an update of the 
BABAR analysis in Ref. and is based on a larger data 
set with improvements in data reconstruction. The data 
were recorded at the T(45') resonance with the BABAR 
detector at the PEP-II storage ring at the Stanford 
Linear Accelerator Center. The integrated luminosity is 
81.9fb^^, corresponding to 89 million BB pairs. 

At the BABAR detector, a five-layer silicon vertex 
tracker (SVT) and a 40-layer drift chamber (DCH), in 
a 1.5-T solenoidal magnetic field, provide detection of 
charged particles and the measurement of their momenta. 
Electrons are detected in a Csl electromagnetic calorime- 
ter (EMC), while muons are identified in the magnetic 
flux return system (IFR), which is instrumented with 
multiple layers of resistive plate chambers. A ring- 
imaging Cherenkov detector (DIRC) with quartz radi- 
ators provides charged-particle identification. 

We fully reconstruct B^ J/ipK^ decays, where 
= TT^ or if^, from the combination of a J/i/; can- 
didate and a charged track K^. The J/ij} candidate is 
reconstructed via a J/?/; e+e" or J/ijj decay 
and is constrained to the nominal J/ijj mass '71] . The elec- 
tron candidates are combined with reconstructed pho- 
tons in the calorimeter to recover some of the energy lost 
through bremmstrahlung. Details of the J/il) reconstruc- 



tion are given in Ref. [2]. Depending on the final state 
of the charmonium meson, the B^ candidates are di- 
vided into two categories, Bee or B^^. The distribution 
in the angle 9i in the J/4> rest frame between one of the 
daughter leptons £ of the J/ip and the line of flight of the 
recoiling is different for signal and background. The 
background peaks for | cos 9i \ near one while the signal 
follows a sin^ 6*^ distribution. We require | cos 6*6 1 < 0.8 
for Bee candidates and | cos6';i| < 0.9 for B^^ candidates. 

Signal yields and charge asymmetries are determined 
by an unbinned maximum likelihood fit to the data. A 
vertex constraint is applied to the reconstructed tracks 
before computing the kinematic quantities of the B* can- 
didate. The beam energy-substituted mass rngs is de- 
fined as 

mES - VGV2 + P-Pb)V^^'-|PbP , (1) 

where -y/s is the total energy of the e^e~ system in 
the T(4S') rest frame, and (S, p) and (iJ^jPe) are the 
four-momenta of the e^e~ system and the reconstructed 
B candidate, both in the laboratory frame. The kine- 
matic variable Ai?7r {^Ek) is defined as the difference 
between the reconstructed energy of the B* candidate 
and the beam energy in the T(4iS') rest frame assum- 
ing ~ ■K^{K^). Signal candidates for B* J/ipir^ 
(B^ J/tpK^) peak in toes at the B* meson mass 
and peak in AE-^ (AEk) at 0. Candidates are required 
to satisfy loose requirements on these variables: lAiS^j < 
120 MeV, \AEk\ < 120 MeV and toes > 5.2GeV/c2. The 
kinematic separation is sufficiently good (see Fig. ^ so 
that no explicit particle identification is required on the 
charged hadron h^, thereby simplifying the analysis. 

The selected sample contains 3801 B^^ and 4053 Bee 
candidates. Figure ^a) shows the toes distribution in 
data fitted to the sum of a Gaussian and an empiri- 
cal phase-space function (Argus function jlfll) describ- 
ing the signal and background components, respectively. 
Figure H^b) shows the AEk distribution for data can- 
didates with Toes > 5.27GeV/c^ fitted to the sum of a 
double Gaussian and a polynomial function, describing 
the dominant B^ J/ipK^ signal and the background 
contribution, respectively. 

The background (bkg) from continuum and generic BB 
decays is characterized using events that are outside the 
signal regions (sidebands of the data sample). Candi- 
dates in the toes sideband are defined by the requirement 
5.20 < Toes < 5.27GeV/c^, where the upper fimit is ap- 
proximately four times the experimental resolution below 
the B mass. Candidates in the AEk and Ai?^ sidebands 
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FIG. 1: (a) The jties distribution for the B candidates in 
data. A fit to the sum of a Gaussian and an empirical thresh- 
old function (dashed curve) is superimposed. The fitted reso- 
lution is approximately 2.5 MeV/c^. (b) The AEk distribu- 
tion for the candidates in data with rriES > 5.27 GeV/ c^. 
A fit to the sum of a double Gaussian and a 3rd order poly- 
nomial function (dashed curve) is superimposed. The fitted 
resolution is approximately 10.5 MeV 



are defined by the requirement 42 < \AEk\ < 120 MeV 
and 42 < lAE'^rl < 120 MeV, where the lower limit is ap- 
proximately four times the AE resolution obtained from 
the fit shown in Fig. ^b). 

We maximize the following extended likelihood func- 
tion: 



M 



L 



where j is the index of the event, i is the index of the 
hypothesis {i = TT,K,bkg), Ni is the yield for each hy- 
pothesis, and M is the total number of events in the 
sample. 

The arguments of the probability density functions 
(PDFs) Pi are a discrete variable a that identifies the cat- 
egory of the B candidate (a = 1 for B^e, a = 2 for i?^^), 
and the kinematic observables (AET^,ph,m-Es)i where ph 
is the momentum in the laboratory frame. We as- 
sume the same PDFs for i?+ and B~ candidates. If we 
define ff'=(A£'^,p/i,mEs) and i^'^'^(A£'^,p,i, toes) as the 



PDFs for Bee and i?^^ candidates, we have 



P, = 



(1 



if a 
if a 



(2) 



where rf^ is the fraction of Bee candidates in a given 
hypothesis. In the following we will drop the superscripts 
ee and /Lt/i when not needed. 

The factor Ci{q) is the fraction of candidates with 



charge q in hypothesis i: 

^^^'^^-\ 1/2(1 + A) ifg 
where Ai is the charge asymmetry: 

Nr - N+ 



+1 

-1, 



(3) 



(4) 



The yields Ni , asymmetries Ai , and fractions r^^ are free 
parameters in the likelihood fit. 

Since the measured variables AE^^ and ph are corre- 
lated, we define a new set of variables: 



D 



AE 



K 



AE^ 



7 



•■K 



Pi 



S = {AEK + AE^)/{D-a), 
n D{D/2~a), 

where 7 is the Lorentz boost from the laboratory frame 
to the T{4:S) rest frame and a = 240 MeV is twice the 
maximum |Ai?^| or |ASif| value for the data sample. 
These variables have the property that {AEt^^D) in the 
pion hypothesis, {AEk,D) in the kaon hypothesis, and 
(1^,11) in the background hypothesis are correlated at 
less than the few percent level. Therefore each Pi can be 
written as a product of one-dimensional PDFs: 



PTr{AE^,Ph,mEs) 
PK{AETr,Ph,mEs) 

Pbkg{AEj,,ph,mEs) 



U{AE^)g^{D)K{mEs) , 

fK{AEK)9K{D)hK{mEs), 
fbkgiS)gbkg(n)hhkgimEs) ■ 



The and fx components are represented by dou- 
ble Gaussians, while /ijr and are described by single 
Gaussians. The parameters of and /ijr are constrained 
to be equal to the parameters of Jk and hx, respectively. 
They are free parameters in the likelihood fit and are ex- 
tracted together with the yields. This strategy reduces 
the systematic error due to possible inaccuracies of the 
Monte Carlo (MC) simulation in describing the AE and 
Toes distributions. 

The (Jtt and gK components are each represented by a 
phenomenological function with seven fixed parameters 
estimated from the MC simulation. They follow an ex- 
ponential shape with Gaussian edges. 

The fbkg component is represented by a linear phe- 
nomenological function with fixed parameters estimated 
from the distribution of S for events in the TOes sideband 
(Fig. Ha)). 

The gbkg component is represented by a phenomenolog- 
ical function with twelve fixed parameters, all estimated 
from the distribution of 11 for events in the toes sideband 
(Fig.Hb)). 

The hf)kg component is represented by the sum of an 
Argus function and a Gaussian function, with fixed pa- 
rameters. The shape parameters are estimated from the 
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FIG. 2: The distribution of (a) E and (b) 11 for events in 
the JTiES sideband in data. The curve corresponds to the 
projection of the best fit. 



distribution of mEs for events in both the I^Ek and AiJ^r 
sidebands. The small number of background events peak- 
ing in the tties signal region is due to candidates recon- 
structed from other B J/il^X decays. From detailed 
MC simulations of inclusive charmonium decays we de- 
termine 40 ± 7 peaking background events in our sample. 

The yields determined with the unbinned maximum 
likelihood fit to the data sample are reported in Tabic 
The correlation coefficient between iV^r and Nk is —0.02. 
The probability to obtain a maximum value of the likeli- 
hood smaller than the observed value is 50%, estimated 
by MC techniques. Figure El shows the distributions of 
Ai?7r for the events in the data, compared with the distri- 
butions obtained by generating events with a parametric 
MC simulation based on the PDFs used in the fit. 

TABLE I: Uncorrected yields Ni, fractions of Bee candidates 
and uncorrected charge asymmetries Ai from the fit to 
the data sample. 



i 




rr(%) 


A. 


n 


242 ± 20 


50.f ±4.f 


0.ff7± 0.084 


K 


4538 ± 70 


46.3 ±0.8 


0.028 ±0.0f5 


bkg 


3074 ± 60 


59.6 ±0.9 


0.0f9± 0.020 



Possible biases in the likelihood estimates were inves- 
tigated by performing the fit on simulated samples of 
known composition and of the same size as the data. 
The samples were generated with parametric MC simu- 
lations based on the PDFs used in the fit. There is no 
evidence of bias in the fitted asymmetries, while a less 
than 1% deviation in the fitted yields from the nominal 
values is present. After correcting the yields for the ob- 
served bias, we obtain = 244 ± 20, Nk = 4548 ± 70, 
and a ratio of branching fractions of (5.37 ± 0.45)% with 
an absolute systematic error of 0.11%. The dominant 




AE (GeV) 

FIG. 3: The Ai5^ distribution in data (points) compared with 
the distribution obtained from a simulated experiment (his- 
togram). The distributions for each simulated component in 
the sample, normalized to the fitted event yields, are also 
displayed. 



source of systematic error is the fixed parameters of the 
PDFs, primarily the PDFs that describe the background. 
Other sources of systematic uncertainty, such as differ- 
ences in the reconstruction efficiencies for J/tpw^ and 
J/ipK^ events and inaccuracies in the description of the 
tails of the AE resolution function, are found to be neg- 
ligible. 

The sample that is used to determine the charge asym- 
metries is defined by imposing as a further requirement 
that the charged track has a polar angle in the range 
[0.41,2.54] radians, includes at least 12 DCH hits, has 
a momentum in the transverse plane pt > 100 MeV/c, 
and points back to the nominal interaction point within 
1.5 cm in the transverse plane and within 3 cm along the 
longitudinal direction. For these tracks the difference 
in tracking efficiency between positively and negatively 
charged tracks - primarily pions - has been studied in 
hadronic events by comparing independently the SVT 
and DCH tracking systems. 

The selected sample contains 3902 J/^h^ and 

3696 B'^ — > J/iph^ candidates. From the likelihood fit 
we obtain the charge asymmetries reported in Table 
The correlation coefficient between Att and Ak is —0.003. 
Using MC techniques we estimate that the probability to 
obtain a fitted asymmetry Ak greater or equal to the one 
observed, in the hypothesis of zero asymmetry, is 6.7%. 

We correct the fitted asymmetries for the small ob- 
served difference in tracking efficiency between positively 
and negatively charged tracks, obtaining A-^ — 0.123 ± 
0.085 and Ak = 0.030 ± 0.015. The uncertainty on the 
corrections contributes 0.004 and 0.005 to the system- 
atic error on A-^ and Aki respectively. The asymmetry 
induced by the different probability of and K~ in- 
teractions in the detector material before the DCH is 
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estimated to be —0.004. This value is conservatively as- 
sumed to be a contribution to the systematic uncertainty. 
The uncertainty in the fixed parameters of the PDFs, 
determined by fits to simulated or non-signal data sets, 
contributes 0.001 to the systematic errors on both 
and Ak- 

Summing in quadrature statistical and systematic er- 
rors, we obtain a 90% C.L. interval of [-0.017, 0.263] for 
A.n- and [0.003,0.057] for Ak- 

In conclusion we measure the ratio of branching frac- 
tions 



[5.37 ± 0.45(stat.) ± 0.11(syst.) 



which is consistent with theoretical expectations and 
with previous measurements. We also determine the 
charge asymmetries 

A-n- = 0.123 ±0.085(stat.)±0.004(syst.), 
Ak = 0.030 ± 0.015(stat.) ± 0.006(syst.) . 

Our results are consistent with previous measurements 
but with significant improvement in the precision. 
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